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1.  Introduction/Background 


Mechanical  analysis  of  microelectromechanical  system  (MEMS)  acceleration  switches  has  been 
successfully  modeled  based  off  of  Newtonian  dynamic  equations  of  motion/  Damping  plays  an 
important  part  in  varying  the  response  of  MEMS  acceleration  switches  and  is  an  integral  part  of 
system  modeling.  In  isolated  environments  filled  with  a  working  fluid  that  dissipates  energy, 
damping  between  parallel  and  sliding  plates  that  are  affected  by  compression  forces  and  internal 
friction  cause  energy  loss  and  is  often  referred  to  as  squeeze-film  and  Couette  damping. 

Reynolds  equations  have  been  applied  to  model  the  squeeze-film  air  damping  in  MEMS 
switches.  The  Reynolds  number  for  MEMS  is  especially  small  due  to  extremely  small 
geometries  and  inertial  forces  are  often  negligible,  so  the  behavior  of  fluids  is  based  primarily  on 
the  viscous  effects. 

Experimental  evaluations  of  MEMS  systems  required  modeling  damping  by  evaluating 
displacement  or  oscillation  using  piezoelectric,  electromagnetic,  or  electrostatic  (capacitive 
sensing)  transduction  mechanisms. A  previous  study  determined  damping  in  oscillating 
cantilever  tuning  forks  by  linearizing  the  Reynolds  equation  at  a  set  gap  distance.^  Other  studies 
measured  the  oscillation  amplitude  at  resonant  frequencies  and  estimated  the  values  from  the 
logarithmic  decrement  of  those  amplitudes.  ’  Experiments  that  did  not  rely  on  resonant 
frequencies  used  capacitance  measurements  or  pressure  sensors  to  relate  changes  in  capacitance 
or  pressure  directly  to  displacement.  The  analytical  methods  and  experimental  methods  for  these 
studies  all  consisted  of  differentiation  of  measured  displacement  to  calculate  damping  forces. 

Previous  reported  studies  were  only  able  to  determine  damping  by  measuring  displacement.  The 
data  for  isolated  MEMS  acceleration  threshold  sensor  or  switches  are  limited  because  switch 
closure  only  quantifies  single  instances  of  voltage  change  (i.e.,  switch  closure).  Two  novel 
methods  are  discussed  as  non-destructive  methods  for  damping  measurement  of  wafer-level 
packaged  MEMS  contact  switches.  The  first  uses  empirical  acceleration  table  results  and 
modeling  a  second-order-differential-equation  fit  to  mimic  the  MEMS  acceleration  switch 
response.  This  method  iteratively  adjusts  the  damping  ratio  until  the  solution  converges. 

The  second  method  uses  harmonic  vibration  empirical  data  to  solve  the  amplification  ratio 
equation  for  both  natural  frequency  and  damping.  The  MEMS  acceleration  switch  is  vibrated  at 
an  amplitude  where  the  switch  begins  to  make  closure  at  its  natural  frequency.  The  frequency 
span  window  values  where  switch  closure  first  occurs  and  finally  opens  are  the  ends  points  for 
solving  a  nonlinear  amplification  ratio  equation.  This  method  is  constrained  to  solve  only 
underdamped  systems;  however,  it  is  more  advantageous  when  quantifying  single  modal 
damping  values  rather  than  impact  excitation  where  many  modes  are  excited. 
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The  two  damping  measurement  teehniques  presented  are  unique  in  their  applieation  and  eover  all 
damping  ratio  values:  underdamped  and  overdamped  systems.  Damping  characterization  will 
benefit  end  users  by  allowing  a  framework  for  modeling  acceleration  switch  response  in  their 
application  and  help  them  correctly  choose  a  closure-acceleration  value  for  the  MEMS  switch. 

1,1  Equations  of  Motion 

A  linear  shock  table  was  used  to  test  the  sensors,  allowing  the  data  acquisition  device  to  plot 
acceleration  and  voltage  change  over  time.  An  open  switch  indicated  little  or  no  voltage  drop, 
but  at  switch  closure,  the  short  circuit  caused  a  change  in  voltage.  In  Fig.  1,  the  voltage  was 
measured  against  time,  and  the  change  in  voltage  corresponded  to  the  g-force  at  the  time  of 
switch  closure.  The  acceleration  and  voltage  change  was  analyzed  with  the  system  modeled  as  a 
mass  with  a  spring  and  damper  in  a  second-order  differential  equation. 


Fig.  1  Example  of  the  data  taken  by  the  data  acquisition  device 


The  impact  table  will  produce  high  amplitude  shocks  and  will  terminate  the  motion  with  shock 
absorbers.  The  resulting  acceleration  forms  a  half-sine  curve,  which  is  applied  as  a  driving 
function  in  the  equation  of  motion.  The  governing  differential  equation  for  the  motion  of  the 
proof  mass  is  given  by  Eq.  1 

mx(t)  +  cx(t)  +  kx(t)  =  Fsin(a)t)  (1) 

or  Eq.  2  in  displacement  form 

o  F 

x(t)  +  2(a>nX(t)  -F  (x)^x{t)  —  — sin(a)t)  (2) 


where  m  is  the  mass  (kg),  c  is  the  damping  coefficient,  k  is  the  spring  constant  (N/m),  F  is  the 
driving  force  (N)  (defined  by  Newton’s  law  Y.F  =  ma),  ^  is  the  damping  coefficient 


c 

2yjkm 


(3) 


and  (Uji  is  the  natural  frequency  (rad/s) 
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-n  =  J^  (4) 

The  general  form  of  the  solution  to  the  system,  assuming  a  harmonie  driving  function  and  an 
underdamped  system,  is  given  by 

x(t)  =  cos(aji^t)  +  C2  sin(a)dt)]  +  i?  sin(ajt  —  cp)  (5) 

where  and  C2  are  displacement  coefficients  of  the  terms  of  the  transient  solution  obtained  from 
the  initial  conditions  (initial  displacement  and  initial  velocity  is  at  zero),  o)  is  the  harmonic 
frequency  (rad/s),  is  the  damping  frequency  (rad/s) 

rOd  =  (6) 

and  R  is  the  coefficient  of  the  steady  state  solution  Eq.  2  derived  by  substituting  the  steady-state 
component  into  the  equation  of  motion 


and  phase  angle,  q),  is 


(p  -  tan  ^ 


Solving  the  equation  for  and  C2  using  the  initial  conditions  results  in  the  real  part  of  the 
solution 


x(t)  =  e  R  sin(q))  cos^oj^t)  +R 


sin((p)-(jL>  cos((p)] 


+  R  sin(a)t  —  q))  (9) 


The  vibration  and  impact  empirical  data  discussed  in  Section  1.3.2  analyzes  two  forms  of  Eq.  9, 
t  «  1  and  t  »  1. 


x(t  «  1)  =  e  R  sin(^)  cosCm^t)  +R 


(^ti)„sin((p)-(Ucos((p))  . 


sin(a)dt)  +  R  sin((jL)t  —  q)) 


x(t  »  1)  = 


7^sin(mt-(p) 


Dividing  the  output  x(t)  by  the  input  leads  to  the  amplification  ratio  for  a  harmonic  second-order 
system 
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1 


(12) 


Xi 


2 


Equation  12  can  be  solved  for  (U^jand  ^  given  that  Xj,  and  the  respeetive  a)’s  are  known. 


1,2  Squeeze-Film  Damping 

Squeeze-fdm  air  damping,  as  shown  in  Fig.  2,  is  created  by  the  compressing  the  gas  between  two 
surfaces  moving  at  a  eertain  deereasing  gap  distanee  and  creating  a  resistant  force  to  the 
movement.  The  damping  pressure  between  two  plates  consists  of  a  viscous  damping  force  and  an 
elastic  damping  force. 


1  1 


Moving  plate 


Squeeze-film  gas 


Gap  distance  Qiq) 


Fixed  plate 


Fig.  2  Squeeze-film  air  damping  between  a  moving 
plate  and  a  fixed  plate 


In  cases  of  high  frequency  oscillations,  the  elastic  force  dominates  the  viscous  force  and  is 
proportional  to  the  displacement  of  the  plate.  The  coefficient  of  the  elastic  damping  force  kg(a) 
as  determined  by  Bao  and  Yang^  is  described  as 


keia) 


(13) 


where  Pa  is  the  ambient  pressure  (Pa),  A  is  the  area  of  the  plate  (m^),  rj  is  the  aspeet  ratio  of  the 
plate  (for  a  square  plate,  rj  =  1),  and  a  is  the  squeeze  number  given  by 


12iJ.(x)f- 

p  /i2 


(14) 


where  p  is  the  coefficient  of  viscosity  of  the  gas  (Pa-s),  and  I  is  the  characteristic  length  of  the 
plate  (m).  As  shown  in  Eqs.  13  and  14,  greater  force  between  the  two  plates  will  be  caused  by  a 
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higher  squeeze  number,  which  can  be  designedly  caused  by  increasing  characteristic  length  or 
decreasing  the  gap  distance.  The  equation  of  motion  is  then  be  described  as 

mx(t)  +  cx(t)  +  (k  +  kg)x(t)  =  F  sin(a)t)  (15) 

Equation  15  can  be  solved  for  the  damping  coefficient  given  the  input  forcing  function  and 
known  input  design  parameters,  m,  k  +  kg,  and  gap  distance  when  the  acceleration  switch 
closes. 

1,3  Analysis  Methods 

Two  methods  were  developed  that  characterize  the  damping  ratio  for  packaged  MEMS 
acceleration  switches.  MEMS  acceleration  switches  only  provide  a  single  point  of  information  in 
an  applied  acceleration  field — ^when  switch  closure  occurs.  Other  MEMS  devices  use  transducers 
like  capacitive,  piezoelectric,  optical,  and  pressure  to  discern  a  continuous  displacement  signal. 
Knowing  displacement  as  a  function  of  time  allows  for  damping  to  be  calculated  via  the  classical 
ring-down  log  decrement  technique.  Additionally,  the  damping  ratio  for  underdamped  systems 
can  be  predicted  by  the  frequency  response  to  a  harmonic  excitation  using  the  -3  dB  rule.  These 
two  methods  are  frequently  used  to  characterize  damping  for  micro -devices.  When  displacement 
is  no  longer  continuously  known,  for  example,  in  MEMS  acceleration  switches,  traditional 
techniques  to  measure  damping  can  no  longer  be  applied.  The  two  methods  discussed  below  are 
ways  to  measure  damping  when  switch  closure  is  all  that  is  known. 

1,3,1  Analysis  Methods  1:  Acceleration  Table  Damping  Predictor 

The  analytic  method  used  the  solution  for  the  equations  of  motion  Eq.  9  to  compute  the  damping 
coefficient  at  the  time  of  specified  voltage  changes.  MEMS  acceleration  switches  developed  for 
60G  fuzing  impact  environments  were  used  to  evaluate  the  non-destructive  method  for  damping 
measurement.  The  Department  of  Defense’s  Joint  Euze  Technology  Program  (JETP)  provided 
the  funding  behind  the  production  of  the  60G  MEMS  impact  switches.  The  JETP  MEMS  impact 
switches  are  omnidirectional,  consisting  of  two  out-plane  contact  (top  and  bottom)  and  one  in¬ 
plane  contact  for  lateral  or  side  acceleration  events.  Eigure  3  illustrate  the  computer-aided  design 
(CAD)  layout  for  the  60G  MEMS  acceleration  switches.  Eor  the  top  and  bottom  contacts  of  the 
JETP  60G  sensors,  the  squeeze-film  air  damping  was  taken  into  account  due  to  the  rectangular 
plate  design. 
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Fig.  3  MEMS  JFTP  60G  impact  switch  layout 

The  numerical  method  approximated  the  damping  by  first  estimating  a  damping  coefficient  value 
and  then  comparing  the  time  value  at  the  determined  damping  to  the  actual  time  of  switch 
closure.  By  adjusting  the  time  value  until  it  was  equivalent  to  the  time  of  switch  closure  and 
modifying  the  damping  coefficient  value  accordingly,  the  damping  coefficient  value  for  the 
function  was  determined  at  the  instance  of  voltage  changes. 

1,3,2  Analysis  Methods  2:  Harmonic  Excitation  Modal  Damping  Predictor 

To  evaluate  the  harmonic  excitation  modal  damping  predictor  (HEMDP)  method,  underdamped 
50-2000G  MEMS  impact  acceleration  switches  were  tested.  These  switches  were  chosen 
because  they  were  known  to  be  underdamped,  which  is  needed  for  the  prescribed  method.  The 
JETP  60G  impact  switch  was  designed  to  have  higher  damping  and  given  that  it’s  a  relatively 
new  device  the  damping  values  were  unknown  and,  hence,  not  used  for  HEMDP  test  setup. 

The  50-2000G  MEMS  impact  switches  used  in  fuzing  applications  were  mounted  to  an 
inductive  shaker  and  vibrated  at  constant  amplitude  for  sine  sweeps  with  frequencies  ranging 
from  100  Hz-10  kHz.  The  switches  are  omnidirectional;  therefore,  switch  closure  is  dependent 
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upon  switch  orientation  during  vibration.  Figure  4  illustrate  the  50G  omnidirectional  acceleration 
switch  with  its  lid  removed. 


Fig.  4  The  50G  MEMS  acceleration  switch  with  its  lid  removed 

Modal  frequency  separation  is  large  enough  that  modal  coupling  effects  can  be  ignored.  The 
MEMS  impaet  switehes  were  evaluated  both  in  the  ±z  out-of-plane  and  in  the  orthogonal  planar 
direction. 

The  modal  damping  values  are  evaluated  by  using  the  harmonie  exeitation  data  where  the  switch 
closes.  This  frequency  is  constrained  by  two  points:  1)  the  frequency  where  the  switch  first 
closes  and  2)  the  frequeney  where  the  switch  opens.  Using  these  two  known  points  and  also  the 
known  designed  gap  distance  the  MEMS  switeh  has  to  travel  during  elosure  allows  one  to  then 
solve  the  nonlinear  amplification  ratio  Eq.  12  for  the  damping  ratio  and  natural  frequency. 


2.  Experimental  Approach 


The  goal  of  the  experimental  approach  was  to  empirically  characterize  MEMS  acceleration 
switches  on  a  linear  impaet  table  and  an  inductive  shaker,  and  use  these  results  to  determine 
damping  values.  During  an  impact  event,  typically  all  modes  are  excited  resulting  in  a 
superposition  of  modal  displacements.  During  harmonie  exeitation,  a  particular  mechanical  mode 
ean  be  excited  generating  a  very  simple  motion.  Damping  values  determined  from  harmonie 
excitation  will  be  very  specific  to  the  mode,  whereas  damping  values  obtained  from  impact  tests 
will  be  generalized. 
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There  were  many  different  parts  to  the  experimental  system  setup  used  to  test  the  sensors,  and 
two  MATLAB  programs  were  written  to  analyze  the  data.  A  linear  shock  table  was  used  to 
produce  a  high  impact  shock.  Four  channels  were  also  used  to  measure  each  sensor  one  at  a  time. 
Six  JFTP  sensors  were  attached  to  the  mounting  plate  and  moves  at  high  speeds  into  a  cushion. 

The  harmonic  excitation  method  used  an  inductive  shaker  with  accelerometer  feedback  control  to 
hold  a  constant  acceleration  value  through  the  sine  sweeps.  Data  acquisition  records  the  time 
versus  switch  voltage  data  and  correlates  into  frequency  versus  switch  voltage.  Switch  closure  is 
indicated  by  a  non-zero  voltage.  MATLAB  post-processing  programs  were  developed  to  sort 
through  each  data  set  and  obtain  the  relevant  information  needed  to  solve  for  damping  and 
natural  frequency. 


2,1  Linear  Shock  Machine 


The  LSM-100  linear  shock  machine  used  pressurized  air  and  a  piston  arm  to  pull  back  the  plate, 
which  compressed  spring  coils  (Fig.  5).  The  piston  arm  would  hit  the  trigger  and  release  the 
plate,  forcing  the  plate  into  a  damping  pad  at  high  speeds.  The  trigger  could  be  moved  at 
different,  measurable  distances  to  vary  the  acceleration  amplitude.  Data  were  collected  and 
displayed  by  a  WinCAT  Data  Acquisition  System  (GFII  Systems,  Inc)  when  the  acceleration  of 
the  plate  surpassed  a  certain  g-force  threshold  specified  by  the  user. 


Fig.  5  Open  LSM-100  linear  shock  machine  (GHI  Systems,  Inc) 
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2,1,1  System  Setup 


The  system  setup  and  WinCAT  software  was  arranged  to  analyze  one  sensor  at  a  time.  The 
power  supply,  channel  input  connection,  and  piston  arm  control  are  shown  in  Fig.  6.  Four 
channels  connected  to  the  data  acquisition  device  were  used  to  collect  data.  Channel  1  was  used 
for  measuring  acceleration  of  the  table  at  impact  and  Channels  2,  3,  and  4  (Fig.7)  measured  the 
voltage  changes  for  each  contact  (top,  in-plane,  bottom).  For  shorter  durations  with  the  harder 
cushion,  data  were  taken  at  a  frequency  of  125  kHz  for  a  duration  of  4.10  ms.  For  longer 
durations  with  the  softer  cushion,  10-bit  data  were  taken  at  a  frequency  of  50  kHz  for  a  time 
duration  of  10.24  ms. 


Fig.  6  Power  supply  (left),  channel  input  connections  (right),  and  piston  arm  control  (bottom) 
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Fig.  7  Three  channel  (left)  and  power  supply  (right)  connections 

2,1,2  Sensors 

Six  JFTP  sensors  (Fig  8)  were  attached  to  a  voltage  divider  circuit  that  was  connected  to  the  data 
acquisition  device.  The  sensors  were  attached  to  a  mounting  plate  in  different  orientations  using 
screws  and  washers. 


Fig.  8  Six  60G  JFTP  sensors  attached  to  mounting  plate 

2,1,3  Impact  Pads 

Different  pads  (Fig  9)  were  used  to  absorb  the  shock  from  the  impact.  Softer  cushions  were  made 
of  rubber  and  the  harder  cushions  were  made  of  a  harder  plastic. 
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Fig.  9  Two  different  shock  absorber  pads:  hard  plastic  (left)  and  soft  rubber  (right) 

2,2  Harmonic  Excitation  Experimental  Setup 

The  harmonic  excitation  method  was  conducted  at  the  US  Army  Research  Laboratory’s  (ARL) 
Weapons  and  Material  Research  Directorate,  Ballistic  Structures  and  Launch  Dynamics 
Laboratory  in  Adelphi,  Maryland.  The  experimental  setup  uses  an  inductive  shaker  controlled  via 
a  feedback  accelerometer.  The  harmonic  excitation  method  for  modal  damping  characterization 
requires  constant  acceleration  amplitude.  The  data  acquisition  system  measures  acceleration  and 
switch  voltage  as  a  function  of  time  and  correlates  these  data  to  the  driving  sine  sweep  frequency 
ranging  from  100  Hz-10  kHz.  The  data  acquisition  software  had  to  be  modified  such  that  switch 
closure  remained  unfiltered  during  sine  sweeps.  Sine  sweep  software  typically  has  a  travelling 
notch  filter  that  removes  frequencies  outside  the  sine  sweep  frequency. 

Figure  10  illustrates  the  inductive  shaker  used  (left)  where  the  MEMS  impact  switch  are 
mounted  on  the  top  of  the  shaker  (right).  The  printed  circuit  board  along  with  the  MEMS  impact 
switches  and  accelerometer  (right)  are  mounted  to  the  top  of  the  shaker;  power,  and  voltage  are 
recorded  by  the  data  acquisition  system.  Eigure  8  shows  the  circuit  board  developed  to  power  the 
voltage  dividing  circuitry  for  each  MEMS  impact  switch.  A  voltage  of  zero  indicates  an  open 
switch,  whereas  a  voltage  of  +5  V  indicates  a  closed  switch. 
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Inductive  shaker  -  Vibration  frequency  sweep 


Fig.  10  Numerical  method;  experimental  data  (top)  and  initial  calculated  model  (bottom) 

2,3  MATLAB  Programs 

MATLAB  programs  where  developed  for  each  testing  method  describe  in  Section  2.  The  first 
described  MATLAB  program  uses  empirical  acceleration  table  results  and  modeling  a  second- 
order-differential-equation  fit  to  mimic  the  MEMS  acceleration  switch  response.  This  method 
iteratively  adjusts  the  damping  ratio  until  the  solution  converges. 

The  second  method  uses  harmonic  vibration  empirical  data  to  solve  the  amplification  ratio 
equation  for  both  natural  frequency  and  damping.  The  MEMS  acceleration  switch  is  vibrated  at 
an  amplitude  where  the  switch  begins  to  make  closure  at  its  natural  frequency.  The  frequency 
span  where  switch  closure  first  occurs  and  finally  opens  are  the  ends  points  for  solving  a 
nonlinear  amplification  ratio  equation.  This  method  is  constrained  to  solve  only  underdamped 
systems;  however,  it  is  more  advantageous  when  quantifying  single  modal  damping  values  rather 
than  impact  excitation  where  many  modes  are  excited. 

2,3,1  MATLAB  Programs  -  Impact  Table  Damping  Model 

MATLAB  programs  were  written  for  both  a  numerical  and  an  analytical  method.  The  numerical 
method  guessed  some  damping  value  and  used  an  ordered  differential  equation  solver  (ode  4  5) 
to  determine  a  time  of  closure  at  the  gap  distance  between  the  proof  mass  and  the  contact.  If  the 
calculated  closure  time  did  not  match  the  closure  time  in  the  data,  the  program  would  estimate  a 
closer  value  of  damping  using  the  relationship 

=  (16) 

where  k  is  a  refinement  factor.  The  analytical  method  used  the  f  zero  function  and  applied  it  to 
the  real  part  of  the  solution  of  the  equation  of  motion  by  setting  the  whole  function  to  zero.  The 
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program  calculated  the  damping  using  the  gap  distanee  between  the  proof  mass  and  the  eontaet 
and  the  experimental  closure  time. 


An  example  of  the  convergence  method  is  portrayed  in  Figs.  1 1  and  12.  In  Fig.  1 1,  a  graph  of 
aetual  experimental  data  for  a  single  test  is  shown  as  well  as  the  model  for  the  ealeulated 
aeeeleration  input.  The  switeh  eloses  the  gap  when  the  displaeement  is  30  pm.  The  ineorrect 
damping  value  in  the  switeh  empirieal  data  shows  that  switch  closure  happens  at  a  later  time,  so 
damping  is  eontinually  inereased  until  the  time  elosure  value  matehes  that  of  the  original.  In 
Fig.  12,  the  eode  eonverged  to  a  final  damping  value  at  2.041  after  42  iterations,  ehanging  the 
damping  by  a  faetor  of  5. 


Fig.  1 1  Numerical  method;  experimental  data  (top)  and  initial  calculated  model  (bottom) 
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Switch  Voltage  [V]  Switch  Voltage  [V] 


Time  [ms] 


Fig.  12  Numerical  method;  experimental  data  (top)  and  final  calculated  model  (bottom) 

2,3,2  MATLAB  Programs  -  Harmonic  Excitation  Modal  Damping  Predictor 

A  well-known  method  to  measure  damping  for  harmonic  excitation  is  by  directly  measuring 
displacement  and  correlating  that  to  the  input  displacement  to  produce  an  amplification  ratio  as  a 
function  of  frequency.  If  input  acceleration  is  held  constant,  then  the  input  displacement 
amplitude  is  proportion  to  inverse  of  frequency  squared.  Therefore  at  higher  frequencies,  less 
input  displacement  is  needed  to  produce  a  constant  acceleration.  Figure  13  illustrates  that  if 
output  displacement  is  known  for  various  excitation  frequencies  then  the  quality  factor  and 
damping  of  the  system  can  be  directly  determined.  Flowever,  in  the  case  of  MEMS  acceleration 
switches,  the  peak  value  of  the  amplification  ratio  is  unknown  because  of  switch  closure. 
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Switch  Voltage  [V]  Switch  Voltage  [V] 


Frequency  Response  to  Harmonic  Excitation 


Fig.  13  The  -3  dB  method  to  quantify  the  quality  factor  and  damping 
of  a  mechanical  system. 


Figure  14  illustrates  that  switeh  closure  only  indicates  the  frequency  at  which  the  switch  closes 
and  then  subsequently  opens  (left).  Thus,  the  only  information  attained  is  the  two  frequency 
values  where  the  switch  first  closes  and  then  opens  and  the  input  displacement,  which  is  input 
acceleration  amplitude  divided  by  the  frequency  squared.  In  the  case  of  Fig.  14  (right),  these 
values  are  0.7395  and  1.094  (non-dimensional  frequency). 


G  1, 


Freq  (Hz) 
Vibration  sweep 


x=1.094 


1  1.5 

Frequency  ratio,  ^ 


Fig.  14  Switch  voltage  as  a  function  of  frequency  (left)  and  the  information  retain  from  an  amplification  ratio 
standpoint  (right) 
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The  values  shown  in  Figs.  14-15  are  used  below  as  an  example  to  ealeulate  the  damping  ratio 
and  natural  frequency  of  MEMS  acceleration  switches  when  all  that  is  known  is  switch  closure 
and  open  frequencies  and  input  acceleration.  Input  displacement  is  determined  by  dividing  the 
constant  input  acceleration  value  by  the  closure  and  opening  frequencies  which  gives  two  values 
that  Eq.  12  can  be  solved.  MATEAB’s  built-in  function  f  solve  can  be  used  to  solve  a  system 
of  nonlinear  equations  (e.g.,  17  and  18). 


x=0.7395  x=1.094 


Frequency  Response  to  Harmonic  Excitation 


Fig.  15  Switch  voltage  as  a  function  of  frequency  (top)  and  the  closure  frequency  window 

plotted  on  top  of  a  frequency  response  graph  that  assume  a  0.1  damping  ratio  (bottom) 
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3.  Results  and  Discussion 


Data  were  obtained  for  MEMS  acceleration  switches  using  the  two  methods  defined  in  the 
Section  2.  Results  from  the  linear  impact  table  and  harmonic  excitations  method  are  organized 
and  discussed.  The  linear  impact  table  method  uses  acceleration  shock  profile  to  induce  MEMS 
acceleration  switch  closure.  This  method  differs  from  the  harmonic  excitation  method  because 
when  the  MEMS  switch  is  put  into  a  transient  shock  environment  multiple  modes  are  excited. 

The  harmonic  excitation  method  uses  an  inductive  shaker  to  excite  single  modes  of  interest.  The 
harmonic  excitation  method  only  works  for  an  underdamped  system,  which  is  a  limiting  factor 
for  this  developed  testing  method.  Both  testing  methods  revealed  nonlinearity  to  damping  as  a 
function  of  input  amplitude.  This  is  due  to  squeeze-film  effects  amplifying  as  a  function  of 
response  velocity  and  gap  closure. 

3,1  Impact  Table  Damping  Results 

Damping  is  expected  to  change  in  a  manner  related  to  the  change  in  velocity  of  the  impact  test. 
This  is  expected  due  to  gap  closure  of  the  MEMS  acceleration  switch.  The  value  of  the  squeeze- 
film  constant,  Eq.  14,  increases  as  gap  distance  decreases  and  input  frequency  increases.  Both  the 
analytical  method  and  the  numerical  analysis  method  were  expected  to  result  in  similarly 
consistent  values  for  the  in-plane  contacts.  Elsing  both  methods,  the  damping  coefficient  for  top 
contact  movement  was  expected  to  be  consistently  higher  than  in-plane  damping. 

The  analytical  method  for  calculating  damping  was  too  sensitive  to  input  variations;  the 
acceleration,  closure  time,  and  frequency  produced  impractical  damping  values  when  applied  to 
the  equation.  Additionally,  manufacturers  closure  gap  tolerances  are  thought  to  be  inconsistent 
and  have  a  fairly  large  standard  deviation.  Static  centrifuge  testing  needs  to  be  conducted  to 
determine  what  the  acceleration  closure  value  is  for  each  orientation — determining  this  will  give 
better  insight  on  how  the  gap  distances  vary. 

The  numerical  method  proved  to  be  a  better  way  of  determining  damping  because  the  analytical 
method  used  many  assumptions  when  solving  the  characteristic  equation.  One  assumption  was 
that  the  system  was  an  underdamped  system,  so  any  cases  of  overdamped  or  critically  damped 
systems  would  be  inaccurately  modeled  using  the  equations  of  motion. 

In  Eigs.  16  and  17,  there  is  a  general  decreasing  acceleration  pattern  over  increasing  closure  time 
for  both  contacts.  There  are  data  points  that  show  constant  peak  acceleration  over  different 
closure  times. 
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Peak  Acceleration  -  [G]  Peak  Acceleration  -  [G] 


[Top  Contact]  Peak  Acceleration  versus  Closure  Time 
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.16  (Top  contact)  peak  acceleration  vs.  closure  time 

[In-plane  Contact]  Peak  Acceleration  versus  Closure  Time 
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Fig.  17  (In-plane  contact)  peak  acceleration  vs.  closure  time 
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Two  or  three  outliers  above  a  damping  value  of  three  were  considered  anomalies  and  were  most 
likely  due  to  the  proof  mass  bouncing  around  or  chattering  inside  the  sensor  and  aliasing  due  to 
the  125-kHz  sampling  frequency  for  data  acquisition,  seemingly  increasing  closure  time. 

Similar  patterns  occur  with  damping  as  with  acceleration.  In  Figs.  18  and  19,  there  is  a  general 
decreasing  damping  pattern  over  increasing  closure  time  for  both  contacts,  but  there  are  patterns 
of  constant  peak  acceleration  over  different  closure  times.  The  top  contact  damping  was  always 
extremely  high  and  had  inflated  damping  values.  There  seemed  to  be  cases  of  overdamping 
because  the  calculated  squeeze-film  factor  was  extremely  large. 

[Top  Contact]  Damping  versus  Closure  Time 
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[In-plane  Contact]  Damping  versus  Closure  Time 
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Fig.  19  (In-plane  contact)  damping  vs.  closure  time 

In  Fig.  20,  the  top  contact  damping  values  are  inflated  due  to  high  squeeze-film  damping.  For 
both  contacts  (Figs.  20  and  21),  there  are  clear  patterns  of  increasing  damping  with  increasing 
peak  acceleration. 
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Damping 


[Top  Contact]  ^  Damping  versus  Peak  Acceleration 
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Fig.  20  (Top  contact)  damping  vs.  peak  acceleration 

[In-plane  Contact]  ^  Damping  versus  Peak  Acceleration 
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Fig.  21  (In-plane  contact)  damping  vs.  peak  acceleration 
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A  relationship  of  increasing  damping  with  increasing  velocity  is  portrayed  in  Figs.  22  and  23. 
The  relationship  between  damping  and  velocity  were  due  to  the  drag  force  in  the  system.  This 
pattern  also  accounts  for  the  relationship  between  damping  and  acceleration  for  the  same  time 
value.  The  values  for  damping  were  inflated  or  inaccurate  due  to  the  sampling  frequency  and 
large  squeeze-fdm  values.  These  factors  severely  impacted  the  calculated  damping  coefficients, 
but  there  were  clear  relationships  shown  in  the  results. 
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Fig.  22  (Top  contact)  damping  vs.  maximum  velocity 


[Top  Contact]  ^  Damping  versus  Maximum  Velocity 
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[In-plane  Contact]  Damping  versus  Maximum  Velocity 


Fig.  23  (In-plane  contact)  damping  vs.  maximum  velocity 

3,2  Harmonic  Excitation  Modal  Damping  Results 

The  damping  results  from  the  harmonie  exeitation  method  are  shown  in  Tables  1  and  2  for  50G 
and  500G  MEMS  acceleration  switches.  When  first  comparing  the  side  (in-plane)  and  out-of¬ 
plane  (top  and  bottom)  damping  values,  the  results  were  unexpected.  However,  upon  further 
review  of  the  contact  surface  area,  the  side  contact  had  roughly  two  times  surface  area  as  the  out- 
of-plane  (top  and  bottom)  contacts.  Squeeze-film  damping  would  be  more  pronounced  in  the  side 
contact  rather  than  the  top  or  bottom  since  it  has  double  the  contact  surface  area. 

Since  the  top  and  bottom  contacts  are  inherently  identical,  the  difference  in  damping  values  are 
attributed  to  the  manufacturing  gap  tolerance  for  the  top  and  bottom  contact  distances.  The  static 
acceleration  switching  threshold  may  deviate  as  much  as  10%  for  the  top  and  bottom  contacts, 
indicating  variation  in  the  gap  distances  as  a  result  of  different  acceleration  threshold  values. 

The  data  in  Tables  land  2  are  for  fairly  low  input  acceleration  amplitude  ranging  from  5-lOG’s. 
Switch  ringing  was  seen  at  input  acceleration  amplitudes  greater  than  20G’s,  where  lateral  switch 
contact  was  sensed.  This  is  believed  to  be  from  contact  happening  at  such  a  great  velocity  the 
switch  rings,  exciting  other  modes  lateral  to  the  input  acceleration  direction.  Larger  input 
acceleration  amplitudes  often  resulted  in  higher  damping  values,  which  are  not  shown  in 
Tables  1  and  2.  This  is  believed  to  be  attributed  to  velocity  dependent  damping  effects  that  were 
also  seen  during  the  impact  table  damping  characterization  tests. 
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Table  1  The  50G  MEMS  acceleration  switch  -  results,  5-lOG  amplitude  excitation 


Mode  #/Shape 

FEA  modal  (Hz) 

Contact 

Freq 

Mode  2  /  Out-plane 

680 

Top 

653.9312 

0.004161 

Mode  2  /  Out-plane 

680 

Bottom 

735.2356 

0.005367 

Mode  3-4  /  In -plane 

798 

Side 

802.6671 

0.044097 

Table  2  The  500G  MEMS  acceleration  switch  -  results,  5-lOG  amplitude  excitation 


Mode  #/Shape 

FEA  modal  (Hz) 

Contact 

Freq 

Damping 

Mode  2  /  Out-plane 

2156 

Top 

2719.009 

0.004860 

Mode  2  /  Out-plane 

2156 

Bottom 

2703.846 

0.001689 

Mode  3-4  /  In-plane 

2482 

Side 

2817.224 

0.009231 

Mode  3-4  /  In-plane 

2482 

Spin 

3216.026 

0.000007 

4.  Conclusions 


For  MEMS  systems  and  systems  characterized  by  low  Reynolds  numbers,  the  damping  force  of 
bodies  moving  in  a  fluid  has  been  found  to  be  linearly  proportional  to  the  speed  of  the  mass.'' 
This  relationship  is  determined  through  analyzing  the  drag  force,  which  is  proportional  to  the 
velocity  and  damping  of  the  system.  In  plotting  the  maximum  velocity  against  the  calculated 
damping  coefficients,  the  results  support  the  correlation  and  show  a  steady  increase  in  damping 
with  increasing  velocity.  This  shows  that  damping  may  be  quantified  as  a  function  of  velocity 
The  pattern  of  decreasing  damping  with  increasing  closure  time  can  therefore  be 
explained  with  an  increasing  drag  force  for  larger  velocity  and  damping  values.  This  pattern 
similarly  accounts  for  increasing  damping  with  increasing  acceleration;  higher  peak  acceleration 
within  the  same  time  duration  causes  a  greater  maximum  velocity. 

Time-closure  values  may  also  have  been  aliased  due  to  the  frequency  at  which  data  were  taken. 
The  proof  mass  bouncing  around  in  the  sensor  was  an  unmeasurable  and  uncontrollable  factor  in 
the  data,  so  time-closure  values  may  have  been  incorrectly  obtained  from  the  data.  The  bouncing 
also  may  have  closed  the  switch  with  other  contacts,  causing  the  MATLAB  code  to  incorrectly 
read  the  proper  closure  time  for  the  correct  contact. 

The  switching  threshold  was  set  to  0.79  V,  which  may  have  affected  the  time-closure  data.  In  the 
event  of  bouncing  or  almost  incomplete  closure,  the  switching  threshold  may  not  have  been 
exceeded  and  the  data  for  time  of  closure  may  have  been  missed  and  resulted  in  subsequent 
closure  being  detected  thus  greatly  increasing  measured  closure  time. 

In  the  future,  a  more  accurate  convergence  algorithm  may  be  required  to  more  efficiently 
converge  on  the  correct  damping  value.  In  addition,  an  increase  in  sampling  frequency  would 
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more  accurately  detect  closure,  which  would  provide  more  accurate  data  to  be  analyzed.  The 
inflated  values  for  the  top  contact  could  be  lessened  if  gap-dependent  squeeze-film  value  model 
was  used.  The  squeeze-film  constant  in  Bao  was  derived  using  a  small  amplitude  of  vibration 
assumption  at  a  set  gap  distance  and  doesn’t  take  into  account  a  reducing  gap  that  goes  to  zero. 

Determining  the  damping  values  of  an  impact  acceleration  switch  helps  categorize  how  the 
switch  functions.  Factors  such  as  mass  and  spring  constants  are  tailorable  and  measurable,  but 
the  majority  of  the  damping  is  caused  by  the  mechanical  switch’s  interaction  with  the  N2 
packaged  gas  environment.  The  novel  method  developed  to  measure  damping  by  measuring 
closure  time  helps  model  designs  to  further  improve  design  to  more  accurately  close  at  the 
expected  acceleration. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


ARL 

CAD 

DUT 

HEMDP 

JFTP 

MEMS 

A 

F 

K 

R 

c 

Cl  ,  C2 

k 

/Ce(o-) 

I 

m 

Pa 

t 

( 

V 

F 

o 


US  Army  Research  Eaboratory 
computer-aided  design 
device  under  test 

harmonic  excitation  modal  damping  predictor 

Joint  Euze  Technology  Program 

microelectromechanical  system 

area  of  the  plate  (m^) 

driving  force  (N) 

refinement  factor 

coefficient  of  the  steady  state  solution  (m) 
damping 

coefficients  of  the  terms  of  the  transient  solution 

spring  constant  (N/m) 

coefficient  of  the  elastic  damping  force 

characteristic  length  of  the  plate  (m) 

mass  (kg) 

ambient  pressure  (Pa) 
time  (s) 

damping  coefficient 

is  the  aspect  ratio  of  the  plate  (for  a  square  plate,  ??  =  1) 
coefficient  of  viscosity  of  the  fluid  (Pa-s) 
squeeze  number 
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harmonic  frequency  (rad/s) 
damping  frequency  (rad/s) 
natural  frequency  (rad/s) 
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